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Abstract 

General formu las for computing the radiation force exerted on arbitrarily oriented and 
arbitrarily shaped nonspherical particles due to scattering, absorption, and emission of 
electromagnetic radiation are derived. For randomly oriented particles with a plane of symmetry, 
the formula for the average radiation force caused by the particle response to external illumination 
reduces to the standard Debye formula derived from the Lorenz-Mie theory, whereas the average 
radiation force caused by emission vanishes. 

Keywords: Radiation force, Electromagnetic scattering and emission; Nonspherical particles 


*Tel.: +1-212-678-5590; fax: +1-212-678-5622. 

E-mail address: crmim@giss.nasa.gov (M. I. Mishchenko). 


1 



J. Quant. Spectrosc. Radial. Transfer , in press (2001) 

It is well known that a particle illuminated by a beam of light experiences a force called 
radiation pressure. This phenomenon is used in laboratories to levitate and size small particles 
[1,2] and affects the spatial distribution of interstellar and circumstellar dust grains [3-5]. In a 
renowned paper published in 1909 [6], Debye showed that the radiation force exerted on a 
spherical particle is directed along the propagation direction of the incident beam and derived a 
simple formula expressing the magnitude of the force in terms of the particle extinction and 
scattering cross sections and the asymmetry parameter. This Note extends Debye’s analysis to 
arbitrarily oriented and arbitrarily shaped nonspherical particles. In addition, it considers the 
radiation force caused by emission of electromagnetic radiation by a nonspherical particle having 
a finite absolute temperature. 

Let an arbitrary particle be illuminated by a time-harmonic parallel electromagnetic wave. If 
the amplitudes of the incident and scattered fields do not change in time, the force due to radiation 
pressure exerted on the particle averaged over the period 2n/coof the incident wave is 

F[ = J^dS (f M (r)) • h (1) 

(Section 6.7 of [7]), where r is the position vector, f M is the Maxwell stress tensor, the 
integration is performed over a closed surface S surrounding the scattering object, and h is the 
unit vector in the direction of the local outward normal to S. Assume for simplicity that the 
scattering object is imbedded in a vacuum. Then the instantaneous value of the Maxwell stress 
tensor is [7] 

f M (r,0 = e 0 E(r,r) ® E(r,r) + c 2 B(r,r) ® B(r,0 - i[E(r,f) -E(r,r) + c 2 B(r,/) • B(r,r)] 7 

= e 0 E(r,O ® E(r,r) + fi 0 H(r,f) ® H(r,r) - |[e 0 E(r,O • E(r,r) + it 0 H ( r >0 • H(r,/)] h 

(2) 

where E is the electric and H the magnetic field, B is the magnetic induction, e 0 is the electric 
permittivity and p. 0 the magnetic permeability of free space, c = l/ >/ e oM^o ' s the speed of light in 
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a vacuum, I is the identity dyad, and ® denotes a dydic product of two vectors. The time 
average of the Maxwell stress tensor is 


(f M (r)) = jRe e 0 E(r) <g> E*(r) + p 0 H(r) ® H*(r) -^(e 0 |E(r)| 2 + p 0 |H(r)| 


2 7 


. (3) 


where * denotes a complex-conjugate value. It is convenient to choose for S a sphere centered at 
the scattering object and having a radius r big enough to be in the far-field zone so that the 
scattered wave is spherical [7]. Then Eq. (1) becomes 

F, = r 2 f dr(f M (rr)Vr, (4) 
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where r = r/r. 

The total electric and magnetic fields at r are vector sums of the respective incident (“inc”) 
and scattered (“sea”) components given by [8] 


E mc (r) = Eg expO'for • r) 
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H inc (r) = — - — V x E inc (r) 
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E sca (r) = — Ef a (r)-r = 0, 
r— >°° r 
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H sca (r) = — - — V x E sca (r) 
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where i = V- I, k = co^/eoPo is the wave number in the surrounding medium, and n wc is the 
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unit vector in the incidence direction. Because the incident and scattered fields are transverse, 
the first and second terms in square brackets on the right-hand side of Eq. (3) do not contribute to 
the integral in Eq. (4). We thus have 

-2 f, 1 2 
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(9) 


The terms 


E mc (r) 


and 


H mc (r) 


are constants independent of r, and their contribution to Ft is 


simply zero. The contribution of the remaining terms follows from Eqs. (5)— (8) and the vector 
identity (a x b) • (c x d) = (a • c)(b • d) - (a • d)(b ■ c) : 
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(10) 


Recalling the extinction theorem [8] and the definitions of the extinction cross section C ext , 

d C . 

differential scattering cross section — ^ca ^ and phase matr j x z [10], we finally derive 
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+ Z l3 (r,n im )U inc + Z I4 (f,n mc )V rinc ] , 


where 7 ,nc = — I— |E n | 2 is the intensity of the incident beam. 

21/Ho 1 ' 

Although the first term on the right-hand side of Eq. (11) represents a force component in 

the direction of n ' nc , the direction of the total force is, in general, different from the direction of 
propagation of the incident beam and depends on its polarization state because of the second 
term. The projection of the total force on any direction n is simply the dot product Fj n . In 
particular, the component of the force in the incidence direction is 
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F, • « inc = — C ext / Inc 


rinc 
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is the asymmetry parameter. 
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is the total scattering cross section, and the quantity 
Cpr = C ext — C sca (cos 0) 


(14) 


(15) 


is the cross section for radiation pressure. 

Although being the result of a lengthy rigorous derivation, Eq. (12) allows a transparent 
physical interpretation. A beam of light carries linear momentum as well as energy. The 
direction of the momentum is that of propagation, while the absolute value of the momentum is 
energy/(speed of light). Since the total momentum of the electromagnetic field and the scattering 
object must be constant, the radiation force exerted on the object is equal to the momentum 
removed from the total electromagnetic field per unit time. Consider the component of the force 
in the direction of incidence. The momentum removed from the incident beam per unit time is 


C ext / inc /c. Of this amount, the part due to absorption is not replaced, whereas the part due to 

scattering is partially replaced by the contribution due to the projection of the moment of the 
scattered light on the direction of incidence. This contribution is equal to the integral of 

/ sca cosO/c over all scattering directions, or / ,nc C sca (cos0)/c . Note that van de Hulst [11] 


used these arguments as a heuristic derivation of Eq. (12). 

If the absolute temperature of the particle T is above zero, then electromagnetic radiation 
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emitted by the particle causes an additional component of the radiation force. The emitted 
radiation is incoherent and does not interact with the incident and scattered radiation, thereby 
generating an independent radiation force component. Furthermore, emission is analogous to 
scattering in that it generates radiation propagating radially in all directions. Therefore, we can 
write the emission component of the radiation force by analogy with the scattering component 
given by the second term on the right-hand side of Eq. (11): 

F 2 (7-) = - — f dcof dr?K el (f,7\co), (16) 

c Jo J 471 

where K el (f,7\co) is the first element of the Stokes emission vector. K el (r,r,0)) is the amount 
of electromagnetic energy emitted by the particle in the direction r per unit solid angle per unit 
frequency interval per unit time and is expressed in terms of the particle extinction, K, and 
phase, Z, matrices as follows [10]: 

K el (r,r,C0) = / b (r,(0)K 11 (f,0))-/ b (7’,0))f dr Z u (r,r,(£>) , (17) 
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(18) 


is the Planck blackbody energy distribution [9]. [Note that the particle is assumed to be 
isothermal, which should be appropriate for micron-sized particles exposed to radiation sources 
of moderate strength.] Unlike F[ , the emission component of the radiation force depends on the 


particle temperature. 

The total radiation force exerted on the particle is the vector sum of the component caused 
by the particle response to the external illumination and the component caused by emission: 
F = F, + F 2 . This force depends on the direction and polarization state of the incident light and 
on the particle orientation, shape, size, refractive index, and temperature. The optical 
characteristics of the particle entering Eqs. (11) and (16) can be computed using one of the 
existing theoretical techniques based on solving the macroscopic Maxwell equations or their 
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integral counterparts [10]. 

Note that for particles in a vacuum, the particle temperature is linked to the incident 
radiation because the total absorption of electromagnetic radiation by a particle must be balanced 
by the total particle emission. In this sense, the forces Fj and F 2 are coupled - both arise from 
the incident radiation. 

Consider now a small volume element comprising randomly oriented particles with a plane 
of symmetry (such as spheres or spheroids) or particles and their mirror counterparts in equal 
numbers and in random orientation. Following the analysis of symmetries performed by van de 
Hulst [11], it is straightforward to show that the emission component of the radiation force 

vanishes, whereas the component caused by the external illumination is now directed along n inc , 
is independent of the polarization state of the incident beam, and is given by 


(F,)=i«‘" c / i "'[(C„ t )-(C sca ){co S 0)] 
where (C ext ), (C sca ), and (c pr ) are the orientation-averaged 


(18) 


extinction, scattering, and 


radiation-pressure cross sections, respectively. This formula is identical to that derived by 
Debye [6] for homogeneous spherical particles. 

Finally we note that the radiation pressure caused by external illumination is accompanied 
by the radiation torque exerted on the particle and given by 


= -J dS r r • [yf M (r)^ 


xr 


(19) 


r 

/ ~ \ 

dr r * 

[(T'M^xrj 
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(cf. page 288 of [7]), where r is the radius of a sphere S centered inside the scattering particle and 


having its surface in the far-field zone. Since r ■ I xr vanishes identically, only the first two 
terms in square brackets on the right-hand side of Eq. (3) contribute to the integrals in Eq. (19). 
The evaluation of this contribution is complicated because it requires the knowledge of not only 
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the transverse component of the scattered electric and magnetic fields, but also of the 
longitudinal component, which we have so far neglected because it decays faster than 1 jr. 

Marston and Crichton [12] computed T, for homogeneous and isotropic spherical particles, 

whereas Draine and Weingartner [3] derived a formula for T[ in the framework of the discrete 
dipole approximation. Emission generates an additional component of the radiation torque 
independent of . 
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